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DESIGN  METHOOOLOGY  FOR  A  LIGHTWEIGHT, 
RESONANT-FREE  PLATEN  FOR  VIBRATION  TESTING 


1.  INTRODUCTION 

Prudent  economics  for  reliability  testing  demands  that  as  large  a 
number  of  test  Items  as  possible  be  placed  In  an  environmental  chamber.  This 
permits  a  rapid  accumulation  of  test  hours,  thereby  gaining  high  statistical 
confidence  for  a  given  test  chamber  time.  However,  serious  difficulties  can 
arise  when  vibration-induced  environmental  stressing  Is  required.  If  frequen¬ 
cies  are  Induced  that  excite  the  natural  frequency  of  portions  of  the  platen  on 
which  the  specimens  are  attached,  a  resonant  condition  results.  Those  test 
specimens  in  the  resonating  area  of  the  platen  experience  a  much  higher  vibra¬ 
tional  stress  than  either  the  test  plan  requires  or  those  items  attached  to  a 
nonresonating  portion  of  the  platen  receive.  In  addition,  the  location  of  the 
resonating  portion  of  the  platen  changes  as  different  frequency  inputs  excite 
different  spring-mass  systems  at  their  individual  natural  frequencies. 

Naturally,  If  at  a  certain  frequency,  one  test  specimen  is  subjected 
to  greater  vibration  stress  than  another,  the  test  may  provide  erroneous  failure 
data  for  test  and  reliability  analysis.  A  problem  such  as  this  can  defeat  the 
entire  purpose  of  the  test  and  cause  confusion  as  to  the  correct  classification 
of  a  vibration-caused  failure. 1 

A  previous  technical  report2  discussed  two  methods  of  solving  or 
minimizing  the  resonance  problem.  One  was  to  design  the  platen  thick  enough 
so  that  the  natural  frequency  of  a  loaded  corner  of  a  platen  would  exceed  the 
highest  frequency  excited  by  the  shaker.  If  this  was  impossible  because  of  a 
weight-displacement  limitation  of  the  shaker,  then  alternatives  in  the  sand¬ 
wiching  of  damping  materials  were  offered  as  an  alternate  solution.  After 
reanalyzing  the  contents  and  causes  of  platen  resonances3  and  beam  deflection 
parameters, 3  it  appears  that  the  design  can  be  refined  to  get  higher  natural 
frequencies  at  highly  significant  reductions  of  weight. 

This  Improvement  can  take  place  by  removing  material  from  outer  por¬ 
tions  of  the  platen,  where  It  Is  a  dead  mass,  and  placing  It  more  in  the  center, 
where  It  adds  strength.  In  reality,  the  platen  takes  the  form  of  an  inverted, 
truncated  pyramid  (ITP)  with  a  square-rectangular-parallelepiped  base  and  top, 
as  Illustrated  in  Figure  1.  This  report  will:  (a)  present  the  equations  and 
methods  for  utilizing  this  type  of  design;  (b)  compare  its  mass  efficiency  with 
that  of  a  constant  thickness  platen;  and  (c)  discuss  the  advantages  of  using 
magnesium  as  the  platen  material.  Appendix  A  contains  the  derivation  of  the 
natural  frequency  equation  used.  Appendix  B  contains  a  computer  program  written 
in  Applesoft  Basic  to  ease  the  burden  of  repetitive  calculations. 

2.  BENDING  MODES  OF  A  RESONATING  PLATEN 

The  Spring-Mass  System. 

When  any  structural  system  vibrates  such  that  the  mass  for  inertial 
force  is  identical  to  the  spring  (of  stiffness)  force,  the  system  vibrates  at 
its  natural  frequency.  On  a  complex  structure  this  can  occur  at  any  number  of 
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Figure  1.  Inverted  Truncated  Pyramid  Type  of  Platen  with 
Square-Rectangular-Parallel  piped  Base  and  Top 
and  Attached  to  Cylindrical  Shaker  Head 


frequencies  when  different  components  combine  or  separate  into  different  subsys¬ 
tems,  each  with  its  own  natural  frequency.  Experience  has  shown,  however,  that 
when  a  test  platen  is  loaded  with  many  items  to  be  vibrated  on  a  hydraulic  or 
electrodynamic  shaker,  the  corners  will  go  into  resonance  starting  at  a  line  of 
bending.  These  lines  of  bending  occur  where  there  is  a  relatively  sudden  and 
significant  increase  in  the  platen’s  area  moment  of  inertia. 

3.  SUDDEN  STIFFNESS  CAUSED  BY  CHANGE  IN  AREA  MOMENT  OF  INERTIA 

Usually  test  platens  are  fabricated  as  a  constant  thickness  slab  of  a 
lightweight  metal,  as  illustrated  In  Figure  2.  There  is  a  rapid  increase  in  the 
section  modulus  where  the  platen  bolts  onto  the  circular  shaker  head,  caused  by 
both  the  additional  thickness  and  higher  modulus  of  elasticity.  Although  the 
edge  of  the  shaker  head  is  not  a  true  fixed  clamp,  in  that  it  fixes  an  entire 
line  of  bending  in  a  vise-like  grip  to  form  a  cantilevered  triangular  beam,  the 
thickness  of  the  platen  forces  a  straight-line  vibration  behavior.  This  has 
been  confirmed  by  experimental  results. 2 

The  ITP  type  of  platen,  as  shown  in  Figure  1,  will  behave  the  same  as 
the  constant  thickness  platen  when  it  experiences  a  sudden  increase  in  area 
moment  of  inertia.  However,  two  important  physical  limitations  must  be  taken 
into  account: 

•  Prudent  fabrication  techniques  require  that  the  base  of  the  platen 
have  a  square  cross  section  instead  of  a  round  cross  section  to  mate  with  the 
shaker  head. 

•  The  shaker  head  is  usually  about  an  inch  thick.  If  the  slopes  of 
the  sides  of  the  inverted  pyramid  were  to  be  projected  to  an  apex  Instead  of 
truncated,  it  would  likely  project  deeper  than  the  thickness  of  the  shaker  head. 

This  indicates  that  the  truncation  of  the  ITP  should  then  have  a  square 
base  with  the  same  thickness  as  the  removed  tip  of  the  pyramid  (the  precise  solid 
geometric  name  for  this  shape  would  be  square-rectangular  parallel  piped).  Also, 
since  the  base's  square  corner  does  not  present  a  sudden  or  rapid  increase  in  the 
area  moment  of  inertia,  the  line  of  bending  is  less  precise.  Therefore,  if  the 
distance  "C“  from  the  platen's  corner  to  the  line  of  bending  is  assumed  to  be  to 
the  shaker-head's  circular  edge,  the  most  conservative  assumption  would  be  used 
In  the  calculations.  In  other  words,  the  actual  natural  frequency  experienced  in 
testing  should  not  exceed  that  used  in  any  calculations. 

Platen  Design  Methodology 

Appendix  A  contains  the  derivation  of  the  equation  for  determining 
the  natural  frequency  of  a  corner  of  an  ITP  type  of  platen.  It  is  Equation  A-27 
on  page  31: 

i - : - 


iE  tan^ 


W 

Cptan$  +  12  - 
A 


where 


ffj  =  Natural  frequency  in  Hertz 

£  =  Young's  Modulus  of  Elasticity  of  Platen  Material, 
lb/sq  in. 

4  *  Angle  of  the  slope  of  the  pyramid  measured  from  the  base 
along  the  edge  formed  by  any  two  sides  that  come  together 
up  to  the  apex  In  degrees  (Figure  1). 

W  =  Weight  of  the  corner-most  test  specimen,  fixture/adapters, 
and  bolt  heads  in  pounds. 

A  *  Platen  area  displaced  by  the  corner-most  test  specimen  and 
its  fixturing  in  square  inches. 

p  =  Density  of  platen  material  in  lb/cu  in. 

C  *  Distance  from  platen  corner  to  shaker  head  in  inches. 

This  equation  will  determine  the  natural  frequency  of  a  platen  corner 
for  various  angles  $.  (Test  fixture  engineers  should  already  have  determined 
all  other  inputs  to  the  equation.)  However,  several  very  practical  consider¬ 
ations  are  not  taken  into  account  by  following  this  equation: 

•  The  edges  around  the  perimeter  of  the  platen  form  a  wedge.  There 
will  probably  be  insufficient  material  to  bolt  the  test  specimens  and  their 
respective  fixtures/adaptors  to  the  platen. 

•  As  was  discussed  previously,  the  shaker  head  may  not  be  of  sufficient 
thickness  to  provide  a  significant  enough  increase  in  moment  of  inertia.  If  it 
is  not,  the  platen  will  have  to  compensate  by  having  its  truncated  portion  take 
on  a  geometric  shape  that  allows  bolting  onto  the  shaker  head. 

«  The  angle  4  was  convenient  for  formula  derivation,  but  should  not 
be  used  on  fabrication  drawings.  The  slope  will  have  to  be  restated  in  terms 
of  the  angle  formed  by  any  one  of  the  pyramid's  sides  with  respect  to  the  base. 
This  conversion  is  simply: 


0  *  ARCTAN  ( 77  tan$)  (2) 

The  problem  of  the  corners  of  the  platen  forming  a  wedge  that  is  too 
thin  can  be  solved  by  adding  a  mass  of  constant  thickness  to  the  top  of  the 
platen.  This  was  not  taken  into  account  during  the  original  derivation,  because 
there  is  a  transition  of  the  neutral  axis  from  the  top  constant  thickness  por¬ 
tion  of  the  platen  to  the  variable  thickness  portion  when  the  line  of  bending 
moves  from  the  corner  towards  the  shaker  head.  The  differential  equation 
becomes  extremely  unwieldy,  and  the  final  equation  would  be  several  magnitudes 
larger  than  Equation  (1).  It  was  decided  that  a  short  iteration  process  would 
perfom  the  job  faster  and  easier,  be  almost  as  accurate  (any  errors  would  be  on 


the  conservative  side,  i.e.,  a  natural  frequency  slightly  higher  than  the  calcu¬ 
lations  indicate),  yet  much  lighter  than  a  platen  of  constant  cross  section. 

The  computer  program  in  Appendix  B  performs  all  of  the  calculations  necessary  to 
design  the  platen.  The  program  performs  in  the  following  order. 

a.  It  steps  through  increasingly  larger  angles  of  4  using 
Equation  (1)  until  the  natural  frequency  is  equal  to,  or  just  exceeds,  the 
maximum  frequency  required  by  the  test  plan  (Figure  3A) . 

b.  The  program  then  adds  a  constant  thickness  mass  to  the  top  of  the 
platen  to  ensure  sufficient  thread  depth  for  mounting  test  specimens.  This 
will  lower  the  natural  frequency  below  that  desired  because  of  the  added  mass 
(Figure  3B). 

c.  Next,  the  program  calculates  a  new  value  for  the  distance  from 
the  corner  of  the  platen  to  the  shaker  head  by  projecting  the  angle  $  to  form  a 
new  wedge  with  the  added  mass  (Figure  3C) .  This  will  drop  the  natural  frequency 
again. 

d.  Using  Equation  (1)  again,  the  program  steps  through  increasingly 
greater  values  of  $  until  the  natural  frequency  is  equal  to  or  slightly  greater 
than  that  of  the  test  plan  (Figure  3D). 

e.  The  program  now  moves  the  distance  from  the  corner  of  the  platen 
back  to  where  it  was  (paragraph  b  above)  with  the  top  constant  thickness  mass, 
but  with  the  value  of  $  as  new*y  calculated  in  paragraph  d.  This  will  raise 
the  natural  frequency  slightly  above  that  required  by  the  test  plan. 

f.  The  program  now  computes  a  new  angle,  0  for  fabrication.  It  is 
the  angle  the  sides  of  the  pyramid  make  with  the  base  (Figure  3E). 

g.  The  minimum  thickness  of  the  mating  surface  of  the  platen  to  the 
shaker  head  is  calculated  by  projecting  the  angle  0  to  the  center  of  the  platen. 
The  distance  from  the  platen's  top  surface  to  this  projected  point  is  the  rec¬ 
ommended  thickness  of  the  platen  at  the  mating  surface  (Figure  3F).  This  will 
assure  that  the  line  of  bending  will  occur  where  the  calculations  assume  it 
will,  at  the  point  where  the  platen's  diagonal  intersects  the  shaker  head. 

h.  Finally,  it  calculates  the  estimated  platen  weight.  If  the  platen 
weight,  when  added  to  the  specimen  and  fixture  weight,  exceeds  that  which  the 
shaker  armature  can  support  for  the  maximum  displacement  required,  then  either 

a  new  layout  with  less  test  specimens  should  be  considered  or  a  platen  material 
of  high  damping  should  be  chosen.  If  the  latter  route  is  taken,  the  calculations 
should  be  rerun  with  a  lower  natural  frequency  until  an  acceptable  weight  is 
reached. 

If  the  computer  program  is  used  for  the  above  set  of  calculations, 
the  display  on  the  computer  screen  will  show  the  S’de  and  bottom  views  of  the 
platen  with  all  critical  dimensions.  At  the  bottom  of  the  screen  is  a  window 
showing  all  of  the  critical  dimensions  required  to  fabricate  the  item,  with  the 
exception  of  the  specimen  mounting  holes.  These  should  have  been  determined 
in  advance  before  the  computer  program  was  run. 


Figure  3A.  Side  View,  45°.  Showing 
the  Calculated  Value  of 
«  for  a  Given  Natural 
Frequency  for  the  Firjt 
Iteration* 


Angle  e  Up  to  the  Plane 
of  the  Top  Surface  to  Font 
a  New  Inverted  Pyraaiid* 


Figure  38.  Side  View,  45°,  Showing 
the  Addition  of  Material 
to  the  Top  Surface  for 
Thread  Depth  for  Bolting 
Test  Specimens  to  Platen* 
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Figure  3D.  Side  View,  45*.  Showing  the 
New  Larger  Calculated  Value 
of  e  for  the  Given  Natural 
Frequency  for  the  Second 
Iteration* 


Figure  3E.  Side  view  Showing  the 

True  Angle  8  Foneed  by  the 
Inverted  Pyrauld'i  Sides 
and  Its  Base 


Figure  3F.  Side  View  Showing  the 

Formation  of  a  Base  to  Mate 
with  the  Shaker  Head  with 
a  Thickness  Equal  to  the 
Projected  Pyrauld  Tip 


NOTE:  Refer  to  Figure  1  for  a  detailed  illustration  and  explanation 
of  the  45°  side  view. 
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SPECIMEN  DAMPING  CAPACITY 


Background . 

Damping  is  the  term  given  for  the  conversion  of  strain  energy  to  heat 
energy  when  the  strain  energy  results  from  the  kinetic  energy  bending  the 
mechanical  system.  Damping  results  from  the  intermolecular  friction  that  is 
experienced  when  there  Is  alternating  compression/tension  strain  from  vibra¬ 
tion.  Without  damping,  every  mechanical  device  that  Is  excited  at  Its  resonant 
frequency  would  fail  from  metal  fatigue.  This  Is  because  the  spring  force  is 
equal  In  magnitude  and  In  the  opposite  direction  from  the  mass  force,  thereby 
putting  the  two  forces  In  dynamic  equilibrium.  The  only  structural  force 
available  to  oppose  the  vibrating  input  forces  (which  are  neither  a  spring 
force  or  a  mass  force,  but  are  external  to  the  mechanical  system)  would  be  the 
damping  force,  which  is  the  product  of  a  damping  constant  multiplied  by  the 
sinusoidal  velocity  of  the  vibration  input. 

Although  the  actual  damping  constant,  or  the  damping  force,  may 
possibly  be  of  interest  in  specific  applications,  the  specific  damping  capacity 
(SDC)  Is  of  greater  Interest  to  the  designer  of  test  fixtures.  When  it  is  not 
possible  to  design  a  platen  that  will  not  experience  resonance  somewhere  within 
the  vibration  spectrum  of  the  test  plan.  It  becomes  necessary  to  select  a 
material  that  has  a  reasonable  amount  of  damping  and  yet  not  compromise  other 
Important  qualities  such  as  density,  modulus  of  elasticity,  or  machlnabll ity. 
Therefore,  when  selecting  material,  there  is  a  compromise  among  relative  damping, 
relative  density,  the  modulus  of  elasticity/density  ratio,  and  relative  machin- 
abil ity. 

5.  SPECIFIC  DAMPING  CAPACITY 

The  specific  damping  capacity  (SDC)  has  units  of  percent  per  cycle 
(l/cyc).  The  value  of  the  SDC  increases  as  the  vibration-induced  stress  in  the 
material  increases.  Therefore,  anytime  an  SDC  number  Is  quoted,  a  stress  level 
is  quoted  also.  Since  the  sole  objective  of  a  vibration  platen  is  to  simply 
hold  bolted  test  specimens  while  in  a  vibration  environment,  you  can  reasonably 
expect  the  lowest  stress  level  values  of  the  SDC  to  be  acceptable. 

Close  examination  of  Figure  4  shows  that  a  rough  estimate  of  the 
impact  of  the  SDC  on  the  vibration  amplitude  can  be  calculated  by  using  an 
Inverse  proportion.  An  example  of  this  would  be  Equation  3: 

Xi  (SDCO  -  X2  (SDC2)  (3) 


where 


Xj  *  The  relative  amplitude  at  resonance  of  a  platen  mode 
of  material  #1 

X2  =  The  relative  amplitude  at  resonance  of  a  platen  made 
of  material  #2 

SDCj  =  Specific  damping  capacity  of  material  #1 
S0C2  ■  Specific  damping  capacity  of  material  #2 


Figure  4.  Influence  of  Specific  Damping  Capacity  on 
Theoretical  Resonance  Amplification  Factor 


If  you  were  to  pick  a  reference  material  from  a  past  test  run  with  a 
logged  vibrational  amplitude  at  resonance,  you  could  easily  substitute  the  SDC 
of  any  other  material  to  determine  its  probable  vibrational  amplitude.  However, 
a  ratio  of  X2  to  Xj,  where  Xj  would  be  the  reference  material,  would  suffice  just 
as  wel 1 . 


6.  MAGNESIUM  AS  A  PLATEN  MATERIAL 

Background. 

When  the  test  program  requires  vibration  frequencies  higher  than  test 
economics  or  physical  constraints  will  permit  a  resonance-free  platen  design, 
then  the  selection  of  fabrication  material  with  a  high  SDC  becomes  important. 

Two  structural  materials  that  offer  SDC's  sufficiently  high  for  a  designer's 
consideration  are  usually  considered,  iron  and  magnesium.  Equation  (1)  illus¬ 
trates  the  importance  of  a  high  elastic  modulus  to  density  ratio  in  any  material 
selection.  The  higher  this  ratio,  the  higher  the  resonant  frequency.  Also, 
since  the  total  weight  suspended  on  the  armature  of  the  electrodynamic  shaker 
greatly  affects  the  system's  displacement  (this  is  most  important  at  very  low 
frequencies),  platen  density  becomes  critical.  In  other  words,  the  lighter  the 
material,  the  higher  the  displacement  of  oscillation. 


Table  1  compares  the  densities  and  the  elastic  moduli  to  density 
ratios  of  two  classes  of  cast  iron  to  that  of  cast  magnesium.  Iron  is  about 
four  times  more  dense  than  magnesium,  and  the  elastic  modulus  to  density  ratio 
can  be  as  low  as  less  than  half  that  of  magnesium.  Both  of  these  qualities 
make  cast  iron  ideal  for  stationary  machinery  or  vehicles  where  traction  and 
vibration  damping  are  important  considerations.  However,  a  platen  for  vibration 
testing  requires  the  opposite  qualities,  and  magnesium  is  the  wiser  choice. 

7.  WHICH  MAGNESIUM  ALLOY  SHOULD  BE  USED? 

The  fixture  designer  can  choose  the  magnesium  for  a  platen  from  the 
alloys  available,  the  fabrication  method  (such  as  sand  cast  or  wrought),  and 
heat  treat  methods.  Extensive  research  into  the  effects  of  these  metal lurgical 
factors  in  relation  to  the  SDC  has  provided  these  significant  findings: 

•  The  SDC  has  an  Inverse  relationship  to  alloy  content;  i.e.,  the 
lower  the  alloy  content,  the  higher  the  SDC. 

•  Cast  magnesium,  which  provides  no  grain  orientation,  has  the  highest 
SDC  of  any  fabrication  method.  Any  working  of  the  magnesium  will  have  adverse 
effects  on  Its  use  as  a  good  vibration  damping  material. 

•  Heat  treating  to  Increase  the  tensile  yield  strength  will  lower  the 

SDC.5 


Figure  5  Illustrates  the  alloying  effect  on  the  SDC  of  magnesium  and 
other  structural  materials.  It  also  shows  how  changing  the  bending  stress 
levels  changes  the  SDC.  Particularly  notice: 

•  How  poor  a  material  aluminum  Is.  Because  of  Its  cost  and  avail¬ 
ability,  aluminum  has  long  been  used  as  a  platen  material,  thereby  overstressing 
test  components  greatly  at  resonance. 

•  How  high  the  SDC  is  for  both  unalloyed  magnesium  and  K1A  magnesium, 
but  the  latter  is  able  to  withstand  much  higher  stress  levels.4 

Figure  6  shows  oscilloscope  traces  of  several  sand-cast  magnesium  and 
aluminum  alloys  that  have  been  allowed  to  decay  in  their  vibration  levels  after 
receiving  the  same  Initial  excitation.  Magnesium  K.1  A  has  a  dramatic  effect  when 
compared  with  the  other  materials.5 

The  K1A  Magnesium  Alloy 

As  may  be  suspected,  K1A  is  a  magnesium  with  only  one  alloy.  Its 
chemical  composition  is  shown  in  Table  2.  The  very  low  amount  of  the  single 
alloy,  zirconium,  is  sufficient  to  bring  the  yield  tensile  stress  up  from  a 
little  over  1,000  psi  for  unalloyed  magnesium  to  a  minimum  of  6,000  psi  for  the 
KlA-type  of  magnesium  (Table  3).  This  yield  stress  level  can  be  expected  to 
more  than  adequately  handle  any  anticipated  streses  from  a  vibration  test.  If, 
however,  there  is  any  reason  to  suspect  the  bending  stress  level  in  the  platen, 
a  simple  stress  calculation  at  the  highest  programmed  "g"  loading  can  be  per¬ 
formed.  If  the  stresses  approach  or  exceed  5,000  psi,  the  computer  program 
found  in  Appendix  B  should  be  rerun  with  a  greater  tread  depth  of  the  top 
portion  of  the  platen. 
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Table  1.  Comparison  of  Density  and  Elastic  Modulus 
for  Iron  and  Magnes1um3.6 


Material 


Oenslt 


lb/ In. 3 


X106  psl 


X106  In. 


Grey  Iron 


0.278 


46.76 


i:|:f 

% 


Malleable  Iron  0.278 


Cast  magnesium 


0.0650  to 
0.0665 


89.93 


100.00  to 
97.74 


Flgur*  5.  Soaping  Capacity  of  Soaa 
Sand  Cast  Alloy* 
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Flgur*  6.  Ojclllograa*  Shoving  Oaaping  Capacity  of  Stvoral  Sand  Cost 
Nagnoilua  and  Aluminum  Allayi—  T«*t  Condition*  Conitpnt 
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Table  2.  Chemical  Composition  of  K1A  Magnesium7 


Element 

Percent 

Zi  rconium 

0.40  -  1.0 

Total,  other  elements 

0.30  Max 

Magnesium 

Remainder 

Table  3. 

Mechanical  Property  Requirements 
for  Separately  Cast  Specimens7 

A1  loy 
and 

temper 

Minimum 

tensile 

strength 

Yield  strength  at  0.2t 
offset  or  at  extension 
under  load  indicated 

Elongation 
in  2  inches 
minimum 

_ 

Percent 

K1A-F 

6,000  0.0029 

14 
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APPENDIX  A 


A  DERIVATION  OF  THE  NATURAL  FREQUENCY  EQUATION 
FOR  VARIABLE  THICKNESS  PLATENS 


1.  INTRODUCTION 

When  it  resonates,  each  corner  of  a  shaker's  inverted,  truncated 
pyramid  behaves  as  a  cantilevered,  variable  thickness,  right  triangular  plate 
as  illustrated  in  Figure  1  of  the  text.  If  loading  is  equal  on  all  four  corners 
of  the  test  specimen,  they  will  resonate  at  the  same  frequency.  If  the  loading 
is  unequal,  the  corners  will  not  only  resonate  at  different  frequencies,  but 
the  cantilevered  line  from  which  the  bending  originates  may  shift  as  well. 
Therefore,  any  generalized  equation  must  be  in  terms  of  (a)  the  distance  from 
the  corner  of  the  platen  to  the  line  where  bending  begins;  (b)  equivalent  mass 
loading;  (c)  area  moment  of  inertia  and;  (d)  Young's  modulus  of  elasticity. 

2.  RALEIGH  METHOD  OF  NATURAL  FREQUENCY  CALCULATION 

The  Raleigh  method  of  natural*  frequency  calculation  has  been  found 
to  produce  values  that  agree  closely  with  experimental  findings  if  accurate 
assumptions  are  made  of  the  physical  characteristics  of  the  vibrating  system. 

The  equation  for  this  method  is: 


I  n  dx 

“  =  -fg  — - -  (A-l 

f  / Y2  dx 

where 

<d  »  Natural  frequency,  radians/sec 

Y  =  Equation  for  the  deflection  of  the  vibrating 
system,  inches 

g  «  Acceleration  due  to  gravity,  which  is  a 
constant  of  386.4  in. /sec2 

The  deflection  is  found  by  solving  the  fourth-order  differential  equation:** 


d4Y  WT 

—  ■*  _ 

dX4  El 


(A-2) 


•Hansen,  H.M.,  and  Citenea,  Paul  F.  Mechanics  of  Vibration,  John  Wiley  and 
Sons,  New  York,  NY.  1952. 

**Popov,  E.P.  Mechanics  of  Materials,  Printice-Hal 1  ,  Inc.,  New  York,  NY.  1952  . 
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where : 


Si 


Wj  *  total  loading  on  the  vibrating  system,  Ib/in. 

£  »  Young's  Modulus  of  Elasticity,  lb/in. 3 
I  ■  Area  moment  of  inertia,  in.* 

The  primary  problem  In  solving  Equation  A-2  Is  determining  the  generalized 
expressions  of  Wj  and  I  for  a  specimen  loaded,  inverted  truncated  pyramid, 
right  triangular,  cantilevered  plate.  The  moment  of  inertia  will  be  calculated 
first. 


Equations  for  Moment  of  Inertia  and  Axis  of  Bendinc 


Following  is  the  equation  for  the  moment  of  inertia  of  a  right- 
triangular  plate  with  an  isosceles  triangular  cross  section  where  the  two 
angles  angles  are  and  the  thickness  changes  "x  tan"  as  shown  in  Figures  A-1A 
and  A-18. 

a.  Moment  of  Inertia. 


(width)(thickness)' 


where: 


Width  -  2x 
Thickness  *  x  tan  $ 


(2x)(x  tan  $) 3 


x4  tan3  <(> 


(A- 3 ) 


To  simplify  for  future  calculation 


Let  D 


tan3  $ 


then  I  =  Dx^ 


(A-4) 


( A- 5 ) 
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Flgur*  A- LA. 

G*OMtry  of  th*  Trltnglt 

Foraad  hy  th*  (nt*rs*ctton 
of  Tao  PliiMS  Of  i  Pyrsald 

Tht  dlstsK*  t  Is  Mssurtd  fro*  th* 
ln**rt*d  pyraald's  comtr  toairds  Its 
cpnttr.  TIM  sngl*  •  Is  Mtsurtd  from 
th*  bts*  to  th*  Him  of  1nt*rs*ct1on. 

K  *|  ^ 

1 

X  «n*_ 1  X  J"( 

Flgur*  A-li.  S*oa*try  of  the  Cross-Section 
of  tht  In»*rt*d  Pyrsald  «t 
OlsUnct  X  from  tht  Corner  and 
Perpendicular  to  Figure  A-1B 


b.  Neutral  Axis  of  Bending.  The  centroid  Is  located  halfway  between 
the  left  and  right  edges  o^  the  platen  (along  the  diagonal)  and  one-third  of  the 
way  down  from  the  top  surface;  therefore. 


N  »  (1/3)  x  tan  <t> 


(A-6) 


Calculation  of  Load 


a.  Platen  Load.  The  cross-sectional  area  varies  as  the  distance  Mx“ 
from  the  corner  increases.  The  change  in  cross-sectional  area  is  one-half  the 
function  of  both  the  width  and  thickness: 


d  area  »  1/2  (d  width)  x  (d  thickness] 
//-  area  *  1/2  f /*(  2dx)(  tan  $  dx) 


//da 


'  dx2 


area  »  tan  $  //< 

At  x  *  0,  all  constants  of  integration 
then  a  *  1/2  x  2tan  $ 


(A-7) 
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This,  then  is  the  cross-sectional  area  at  any  distance  "x"  from  the  platen 
corner.  To  get  the  distributed  platen  load,  simply  multiply  the  area  by  the 
density  p  and  obtain: 

Wp  =  1/2  p  x  2tan  <f>  ( A-8 ) 


b.  Specimen  Load.  The  natural  frequency  of  a  cantilever  is  most 
affected  by  masses  that  are  furthest  from  the  edge  of  bending.  The  greater 
the  mass  and  the  further  it  is  from  this  edge,  the  lower  the  natural  frequency. 
This  is  because  the  natural  frequency  is  related  to  the  deflection;  therefore, 
to  simplify  the  deflection  equation,  one  can  assume  that  the  weight  per  displaced 
unit  area  of  the  corner-most  test  item  is  the  same  for  the  entire  triangular 
plate  being  considered.  (Test  prudence  would  dictate,  wherever  possible,  that 
the  lightest  items  per  square  inch  go  at  the  extreme  corners  of  the  platen  to 
get  the  highest  natural  frequency).  Include  the  weight  of  the  test  item,  the 
adapters  (if  any),  and  the  bolt  heads  to  obtain  the  most  accurate  estimate  of 
the  cantilever's  loading.  This  total  specimen  weight  shall  be  designated  as 
"Wx."  Next,  calculate  the  area  displaced  on  the  platen  by  the  test  specimen 
and  all  associated  hardware  and  designate  this  area  as  "A.”  This  provides  a 
load  pressure  of  W/A  in  pounds  per  square  inch  for  the  corner-most  position  on 
the  test  platen.  What  is  now  needed  is  the  equivalent  distributed  load  as  the 
right  triangular  area  expands  from  the  corner  to  the  line  of  bending.  This 
load  increases  at  a  rate  of  "2x"  since  the  width  expands  at  twice  the  rate  as 
the  diagonal  distance  from  the  platen  corner  increases.  The  specimen  load 
becomes : 


W 

Wp  =  2 -  x  ( A-9 ) 

A 


c.  Total  Platen  Load.  The  total  platen  load  is  the  addition  of  the 
specimen  load.  Equation  A-9,  to  the  platen  load.  Equation  A-8.  It  now  becomes: 


Wf  a  Wx  +  Wp 


WT  =  2 - x  +  1/2  p  x2tan  <f>  (A-10) 

A 


To  simplify  for  future  number  manipulation: 


W 

let  B  =  2 —  ( A- 1 1 ) 

A 
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and  G  3  1/2  p  tan  $ 


(A- 12) 
( A- 13) 


then  Wt  =  Bx  +  Gx^ 


With  both  the  load  and  moment  of  inertia  equations  developed,  the  deflection 
can  now  be  derived. 


Calculation  of  Deflection 


From  Equation  A-2  we  get  the  equation  for  deflection: 

d4Y  WT 

dx^  El 

The  solution  to  this  fourth  order  differential  equation,  fortunately,  is  four 
easy  steps  of  integration: 


S  *  J * Wf  dx  +  S0 

( A— 14 ) 

M  =  J*Sdx  +  Mq 

( A- 15 ) 

f  M 

¥  =  /  -  dx  +  ¥0 

( A- 1 6 ) 

J  El 

where: 


■/’ 


Ydx  +  Yf 


(A- 17) 


S  3  shear  load ,  lb 
M  =  moment  load,  in. -lb 
¥  =  slope,  rad 
Y  3  deflection,  in. 


The  mathematical  work  will  follow  these  steps  with  integration  starting  from 
the  free  end  and  going  to  the  fixed  end,  as  shown  in  Figure  A-l. 
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/A? 


Shear  Calculation 


From  Equation  A-14 


From  Equation  A- 13 


S 


/wT  dx  ♦  Sq 


then 


therefore 


WT  *  Bx  +  Gx2 

S  =y*(Bx  -  Gx2)dx  +  SQ 
S  =•  1/2  Bx2  +  ( 1/3 ) GX3  +  S0 
at  X  *  0,  Sq  *  0 


S  =  1/2Bx2  +  (i/3)Gx3 


Moment  Calculation 


(A- 18) 


From  Equation  A- 1 5 


Substitute  the  value  of  S  from  Equation  A-18,  into  the  above  equation: 


M  *  j { l/2Bx2  +  (l/3)GX3)dx  +  M0 
M  *  ( l/6)Bx3  +  (1/12)GX4  +  M0 


at  X  =>  0,  M0  =  0 


Therefore : 


M  =  ( 1/6) Bx3  +  ( 1/12) GX4 


( A- 1 9 ) 
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Slope  Calculation 


From  Equation  A-16 


f- 

»  El 


dx  +  f0 

From  substituting  Equations  A-5  and  A-19  into  Equation  A-16,  we  get: 

1/6Bx3  +  i/12Gx4 


dx  +  f0 


ED  xfl 


B  G 

-  In  x  +  -  x  +  fr 


At  X 
Then 


6  ED 

C  (see  Figures  2  and  3),  f  *  0 


12  ED 


•(- 


B  GC 

6  ED  12  ED 


B  G  /  B 

f  »  -  1  n  x  + -  x  -  I  1  n  C  + 

6  ED  12  ED  \  6  ED 

Deflection  Calculation 

From  Equation  A- 1 7 : 

Y  *  Y  dx  +  Y0 

Substitute  in  the  value  from  Equation  A-20  and  we  get 


GC 


12  ED. 


ff  B 

/  /  _ —  |n  v  x 

G 

B 

— —  1  n  C  - 

GC  \ 

1  1  111  A  ' 

}  \  6E0 

12ED 

6ED 

12ED  / 

(A-20) 


dx  +  Yr 
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xB 

Y  3 

In  x  - 

B 

G 

„  B  1  n 

X2 _ 

C 

-I-  v  _  . 

GC 

_ x  +  y 

6ED 

6ED 

24ED 

6ED 

12ED  ° 

Pij  x  *  C,  Y  - 

0 

then 

CB 

Y0  * - In  C  + 

6ED 

BC 

GC2 

BC 

——In  C  + 
6ED 

GC2 

6E0 

24ED 

I2ED 

BC 

Y0  *  - 

6ED 

GC2 

+ 

24ED 

BC 

Y  *  - 

In  x  - 

B 

G  2 
x  +  — —  x** 

B  In  C 

GC 

8C  GC2 

X  +  — •  +  —  ■  -  — 

°  6ED 

6ED 

24ED 

6ED 

12ED 

6ED  24ED 

1 

Y  a  — 

CO 

X 

In  x 

-  4Bx  +  Gx2 

-  4Bx  In  C 

-  r* 

+  4BC  +  GC2^  (A. 

24ED 

c  O  wA 

This  completes  the  derivation  of  the  deflection  equation  for  a  right  triangular 
plate  with  an  Isosceles  triangular  cross  section. 

Calculation  of  Natural  Frequency 

Equation  A-l  provides  the  Raleigh  equation  for  determining  the  natural  frequency 
of  a  vibrating  system.  It  is: 


hi 


/Ydx 

/Y2dx 


This  requires  an  additional  integration  step  of  the  numerator  and  an  additional 
integration  step  after  squaring  in  the  denominator.  To  prevent  confusion  in 
the  calculations,  let's  set: 


k  =  /Ydx  (A-22) 

and 

j  =  /Y2dx  (A-23) 
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First,  let's  solve  for  "It"  by  substituting  Equation  A-21  into  Equation  A-22: 


*  J_  [(  4  Bx  In  x  -  4Bx  +  Gx2  -  4Bx  In  C-2GCx  +  4BC  +  GC2  ) 

o a rn  I\  / 


24E0  Jx 
0 


dx 


-2Bx2  ♦  1/3Gx3  -  2Bx2  In  C 


k  = 


2Bx2  In  x  -  2Bx2  +  1/3Gx3  .  2Bx2  In  C  -  GCx2  +  4BCx  +  GC2x 


k-  o 


1 

24ED 

k  «  — —  ( 28C2  In  C  -3BC3  *  1/3GC3  -  2BC2ln  C  -  GC3  +  4BC2  +  GC3 
24ED  V 


k  *  — L  (  1/3QC3  +  BC2  ) 

24ED  \  J 

c2  /  \ 

K  *  -  1/3  GC  +  B 

24ED  \  / 

Now,  solve  for  "j"  by  substituting  Equation  A-21  into  A-23: 
/  1 


( A— 2  4 ) 


24E0 


2  /  4  x  B  In  x  -  4Bx  +  Gx2  -  4Bxln  C  -  2GCx  ♦  4BC  +  GC2 )2dx 


j 


- -)2  f  ( 16  B2x2 (In  x)2  -  32B2x21  n 

24ED  /  J  V 


+  8BGx3l n 


32B2x2  In  C 


lnx  -  16BCGx2  In  x  +  32B2  C  x  In  x  +  8BC2Gx  In  x  +  16B2x2  -  8BGx3 
+  32B2x2  In  C  ♦  24BCGx2  -  32B2Cx  -  24BC2Gx  ♦  G2x4  -  8BGx3  In  C 
-  4CG2x3  +  6C2G2x2  +  16B2x2( 1 n  C)2  +  16BCGx2  In  C  -  32B2Cx  In  C 
-8BC2Gxl n  C  -  4G2C3x  +  16B2C2  +  8  BC3G  +  G2C4)dx 
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By  substituting  Equation  A-24  Into  A-22  and  A-25  Into  A-23  and  both  of  these 
Into  Equation  A-l,  Me  get: 


The  natural  frequency  is  usually  expressed  as: 


1 
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Substituting  Equations  A-4 ,  A-ll,  and  A- 12  into  Equation  A-26,  we  get: 


This,  then,  is  the  final  equation  for  determining  the  natural  frequency  of  a 
corner  of  a  truncated  inverted  pyramid.  This  equation  is  found  in  the  main 
text  of  this  report  as  Equation  1. 
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APPENDIX  B 


COMPUTER  PROGRAM  FOR  DETERMINATION 
OF  OPTIMUM  DESIGN  PARAMETERS 


This  program  determines  the  optimum  design  parameters  for  an  inverted 
truncated  pyramid  type  of  platen  for  holding  test  specimens  during  vibration 
testing. 


t.J.  LIPP 
USAAMCCOM 

Product  Assurance  Directorate 
For  Chemical  Systems 
1985 

THE  PROGRAM  REQUIRES  THAT  THE  FIXTURE  DESIGNER  COMPLETE  A  LArO'JT  OF  THE 
TOP  OF  A  SQUARE  VIBRATION  PLATEN  BEFORE  GETTING  ON  THE  COMPUTER. 

THE  DESIGNER  IS  EXPECTED  TO  HAVE: 

1.  ARRANGED  ALL  Of  THE  TEST  SPECIMENS  SO  THAT  THOSE  WITH 
THE  LEAST  HEIGHT  (INCLUDING  FIXTURE  ADAPTERS  AND  BOLTS) 

PER  DISPLACED  PLATEN  AREA  ARE  AT  THE  FOUR  CORNERS; 

2.  THE  TEST  SPECIMENS  AS  CLOSE  TOGETHER  AS  POSSIBLE. 

YET  HITH  ROOM  FOR  HRENCHES,  POWER  AND  SIGNAL  CABLES. 

ETC.  TO  KEEP  THE  PLATEN'S  SQUARE  DIMENSIONS  AS  SHORT 
AS  POSSIBLE; 

3.  SELECTED  THE  PLATEN  MATERIAL  AND  TO  KNOW  ITS  YOUNG'S 
MODULUS  OF  ELASTICITY  AND  DENSITY; 

4.  DETERMINED  THE  MINIMUM  DEPTH  REQUIRED  FOR  THE  BOLT 
THREADS  TO  FASTEN  THE  TEST  SPECIMENS  AND  FIXTURES  DOWN 
TO  THE  PLATEN. 

HAVE  YOU  COMPLETED  THESE  PRELIMINARY  DESIGN  DETAILS?  <Y  OR  N> 

IMPORTANT  NOTE:  ALL  RESPONSES  TO  THE  FOLLOWING  QUESTIONS  ARE  TO  BE 
IN  DECIMAL  FORM,  NOT  FRACTIONS. 


»vw 


what  is  the  length  of  a  side  of  your 

SQUARE  PLATEN  IN  INCHES? 

?42 

WHAT  IS  THE  DIAMETER  OF  THE  SHAKER  HEAD 
IN  INCHES? 

?  1  8 

WHAT  IS  THE  MINIMUM  MATERIAL  THICKNESS 
REQUIRED  TO  BE  TAPPED  TO  BOLT  THE  TEST 
SPECIMENS  TO  THE  PLATEN  IN  INCHES? 

?  .  5 

WHAT  IS  THE  MAXIMUM  FREQUENCY  REQUIRED 
(N  THE  TEST  PROGRAM  IN  HERTZ? 

?  500 

WHAT  MATERIAL  HAVE  YOU  CHOSEN  FOR  THE 
PLATEN? 

7MAGNESIUM  K I  A 

WHAT  IS  MAGNESIUM  KIA'S  TOUNG'  S  MODULUS 
OF  ELASTICITY  IN  POUNDS  PER  SQUARE  INCH? 
77000000 

WHAT  IS  MAGNESIUM  KIA'S  DENSITY  IN 
POUNOS  PER  CUBIC  INCH? 

?0.07 

ASSIGN  EACH  CORNER  OF  THE  PLATEN  A 
NUMBER  SO  THE  WEIGHTS  OF  THE  AREAS  DIS¬ 
PLACED  BY  THE  TEST  SPECIMENS  AND  THEIR 
FIXTURING  CAN  BE  IDENTIFIED  BY  THEIR 
LOCATION  ASSIGNMENT. 


WHAT  IS  THE  WEIGHT  OF  THE  TEST 
SPECIMEN,  FIXTURING,  BOLTS,  WASHERS, 
ETC.  IN  CORNER  I? 

?1  1 

WHAT  IS  THE  AREA  DISPLACED  ON  THE 
PLATEN  BY  THE  TEST  SPECIMEN,  FIXTUR¬ 
ING,  ETC.  IN  CORNER  1 ? 

?!  I 

THE  CORNER'S  LOAD  IS  1  P S I 
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WHAT  IS  THE  WEIGHT  OF  THE  TEST 
SPECIMEN,  FIXTURING,  BOLTS,  WASHERS, 

ETC.  IN  CORNER  2? 

?  1  1 

WHAT  IS  THE  AREA  DISPLACED  ON  THE 
PLATEN  BY  THE  TEST  SPECIMEN,  FIXTUR¬ 
ING,  ETC.  IN  CORNER  2? 

?  1  1 

THE  CORNER' S  LOAD  IS  1  PS  I 

WHAT  IS  THE  WEIGHT  OF  THE  TEST 
SPECIMEN,  FIXTURING,  BOLTS,  WASHERS, 

ETC.  IN  CORNER  3? 

?1  1 

WHAT  IS  THE  AREA  DISPLACED  ON  THE 
PLATEN  BY  THE  TEST  SPECIMEN,  FIXTUR¬ 
ING,  ETC.  IN  CORNER  3? 

?11 

THE  CORNER'S  LOAD  IS  1  PS  I 

WHAT  IS  THE  WEIGHT  OF  THE  TEST 
SPECIMEN,  FIXTURING,  BOLTS,  WASHERS, 

ETC.  IN  CORNER  4? 

?  1  1 

WHAT  IS  THE  AREA  DISPLACED  ON  THE 
PLATEN  BY  THE  TEST  SPECIMEN,  FIXTUR¬ 
ING,  ETC.  IN  CORNER  4? 

?1  1 

THE  CORNER'S  LOAD  IS  1  PS  I 
CORNER  NUM8ER  1  HAS  THE  GREATEST 
LOAD  OF  1  POUNDS  PER  SQUARE  INCH. 

PRESS  <ANY  KEY >  TO  CONTINUE. 

LENGTHS:  L-42  IN.  D-18  IN. 

TO  CONTINUE  PRINTING  DESIGN  INFORMATION 
PRESS  <ANY  KEY>. 

TO  CONTINUE  PRINTING  DESIGN  INFORMATION 
PRESS  <AN Y  KEY > . 

T-.5  IN.  P-10.39  IN.  B-7.79  IN. 

TO  CONTINUE  PRINTING  DESIGN  INFORMATION 
PRESS  <AN Y  KEY>. 
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to! 

V«*2 

* V? 


ANGLE  THETA-40  DEG.  PLATEN  WT.-810  LBS. 
WHEN  FINISHED, PRESS  < ANY  KEY  > . 

DO  YOU  WISH  TO  TRY  ANOTHER  DESIGN? 

IF  SO,  PRESS  KEY  <C>;  IF  NOT,  PRESS 
ANY  KEY. 

GOODBY 

1PRIO 


',■.*** 

*S’Jl 

U 

•  »J  * 

,sj5 


r 


.y\ 


m 

fry 

vyi 
»♦! 
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1  IF  PEEK  (104)  <  >64  THEN  POKE 

103,1 :  POKE  104,64:  POKE  163 
84,0:  PRINT  CHRJ  (4>"RUN  PL 
ATEN" 

5  SPEED*  150 

10  REM  THIS  PROGRAM  DETERMINES 
THE  OPTIMUM  DESIGN  PARAMETER 
S  FOR  A  TEST  PLATEN 
20  REM  FOR  VIBRATION  TESTING  OF 
TEST  SPECIMENS  WITHOUT  REAC 
HING  A  RESONANT 

30  REM  CONDITION  WITH  A  MINIMUM 
DEAD  WEIGHT  LOAD  ON  THE  S 
HAKERS  ARMATURE. 

50  HOME 
60  VTAB  (6) 

70  PRINT  "  THIS  PROGRAM  DET 

ERMINES  THE" 

80  PRINT  "  OPTIMUM  DESIGN  PA 

RAMETERS  FOR" 

90  PRINT  "  AN  INVERTED  TRUNC 

ATED  PYRAMID" 

100  PRINT  "  TYPE  OF  PLATEN 

FOR  HOLDING" 

110  PRINT  "  TEST  SPECIMENS  D 

URING  VIBRAT-" 

115  PRINT  "  ION  TESTING" 

117  PRINT 

120  PRINT  TAB (  14); "BY  L.J.LIPP 

m 

125  PRINT 

130  PRINT  TAB (  1 4 ) ; "U SAAMCCOM" 

140  PRINT  "  PRODUCT  ASSURANC 

E  DIRECTORATE" 

150  PRINT  "  FOR  CHEMICA 

L  SYSTEMS" 

160  PRINT  TAB (  17);"1985" 

170  FOR  X  ■  10  TO  4000:  NEXT  X 
200  HOME 

210  VTAB  (5) 

220  PRINT  "THIS  PROGRAM  REQUIRES 
THAT  THE  FIXTURE" 

230  PRINT  "OESIGNER  COMPLETE  A  L 
AYOUT  OF  THE  TOP" 

240  PRINT  "OF  A  SQUARE  VIBRATION 
PLATEN  BEFORE" 

250  PRINT  "GETTING  ON  THE  COMPUT 
ER.  THE  DESIGNER" 

260  PRINT  "IS  EXPECTED  TO  HAVE:" 


270  PRINT 

280  PRINT  TAB (  5);"1.  ARRANGED 
ALL  OF  THE  TEST  SPECI-" 

290  PRINT  TAB (  8);"MENS  SO  THAT 
THOSE  WITH  THE" 

300  PRINT  TAB (  8);"LEAST  WEIGHT 
(INCLUDING  FIXTURE" 

310  PRINT  TAB (  8);"ADAPTERS  AND 
BOLTS)  PER  OISPLAC-" 

320  PRINT  TAB (  8 ) ; "ED  PLATEN  AR 
EA  ARE  AT  THE  FOUR" 

324  PRINT  TAB (  8 ); "CORNERS ;  " 

330  PRINT 

340  PRINT  TAB (  5);"2.  THE  TEST 
SPECIMENS  ARE  AS  CLOSE" 

350  PRINT  TAB (  8) ; "TOGETHER  AS 
POSSIBLE,  YET  LEAVE" 

360  PRINT  TAB (  8);"R00M  FOR  WRE 
NCHES,  POWER  AND" 

370  PRINT  TAB (  8);"SIGNAL  CABLE 
S,  ETC.  TO  KEEP  THE" 

380  PRINT  TAB (  8);"PLATEN'S  SQU 
ARE  DIMENSIONS  AS" 

390  PRINT  TAB (  8);"SH0RT  AS  POS 
SIBLE;" 

400  PRINT 

410  PRINT  TAB (  5>;»3.  SELECTED 
THE  PLATEN  MATERIAL  AND" 

420  PRINT  TAB (  8); "KNOWS  ITS  YO 
UNGS  MODULUS  OF" 

430  PRINT  TAB (  8 ) ; "EL ASTI C I TY  A 
ND  DENSITY;" 

440  PRINT 

450  PRINT  TAB (  5);"4.  DETERMINE 
0  THE  MINIMUM  DEPTH" 

460  PRINT  TAB (  8); "REQUIRED  FOR 
THE  BOLT  THREADS" 

470  PRINT  TAB (  8);"T0  FASTEN  TH 
E  TEST  SPECIMENS" 

480  PRINT  T AB (  8>;"AND  FIXTURES 
DOWN  TO  THE  PLATEN." 

490  PRINT 

500  PRINT  "HAVE  YOU  COMPLETED  TH 
ESE  PREL I M I  NARY" 

510  PRINT  "DESIGN  DETAILS?  <Y  OR 
N  >  " 

515  GET  AS 

520  IF  AS  <  >  "Y"  AND  AS  <  >  " 

y"  ANO  AS  <  >  "N"  AND  AS  < 

>  "n"  THEN  2000 

530  IF  AS  *  "N"  OR  AS  »  "n"  THEN 
21  00 
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550  PRINT  "IMPORTANT  NOTE:  ALL  R 
ESPONSES  TO  THE" 

560  PRINT  "FOLLOWING  QUESTIONS  A 
RE  TO  BE  IN" 

570  PRINT  "DECIMAL  FORM,  NOT  FRA 
CTIONS." 

580  PRINT 

590  PRINT  "WHAT  IS  THE  LENGTH  OF 
A  SIDE  OF  YOUR" 

600  PRINT  "SQUARE  PLATEN  IN  INCH 
ES?":  INPUT  L 
610  PRINT 

620  PRINT  "WHAT  IS  THE  DIAMETER 
OF  THE  SHAKER  HEAD" 

630  PRINT  "IN  INCHES?":  INPUT  Dl 
A 

640  PRINT 

650  PRINT  "WHAT  IS  THE  MINIMUM  M 
ATERIAL  THICKNESS" 

660  PRINT  "REQUIRED  TO  BE  TAPPED 
TO  BOLT  THE  TEST" 

670  PRINT  "SPECIMENS  TO  THE  PLAT 
EN  IN  INCHES?" 

680  INPUT  TK 

681  PRINT 

682  PRINT  "WHAT  IS  THE  MAXIMUM  F 
REQUENCY  REQUIRED" 

685  PRINT  "IN  THE  TEST  PROGRAM  I 
N  HERTZ?" 

687  INPUT  MF 
690  PRINT 

700  PRINT  "WHAT  MATERIAL  HAVE  YO 
U  CHOSEN  FOR  THE" 

710  PRINT  "PLATEN?":  INPUT  MS 
720  PRINT 

730  PRINT  "WHAT  IS  ";M$;"*S  YOUN 
GS  MODULUS" 

740  PRINT  "OF  ELASTICITY  IN  POUN 
DS  PER  SQUARE  INCH?" 

750  INPUT  E 
760  PRINT 

770  PRINT  "WHAT  IS  ";M$;"'S  DENS 
ITY  IN" 

780  PRINT  "POUNDS  PER  CUBIC  INCH 
?":  INPUT  DEN 

784  REM  THE  FIRST  SUBROUTINE  SE 
LECTS  THE  PLATEN  CORNER  WITH 
THE  HIGHEST  SPECIMEN  LOAD. 
787  GOSUB  2500 

800  REM  DETERMINE  THE  ANGLE  PHI 
FOR  THE  ANGLE  OF  THE  PVRAMI 


D'S  EDGES. 

805  SPEED*  255 
81  0  GOSUB  3000 

820  REM  THE  THIRD  SUBROUTINE  AD 
DS  THE  MATERIAL  THICKNESS  TO 
THE  PLATEN, 

830  REM  AND  CALCULATES  A  NEW  TE 
MPORARY  PLATEN  LENGTH  SO  THE 
EDGES  WILL  BE  WEDGE 
840  REM  SHAPED  AGAIN. 

850  GOSUB  4000 

860  REM  THE  4TH  SUBROUTINE  RECA 
LCULATES  THE  ANGLE  AGAIN,  BU 
T  CALLING  IT  PH  I  2 
870  REM  WHICH  IS  THE  FINAL  ANGL 
E  OF  THE  PYRAMIDS  EDGES. 

880  GOSUB  5000 

890  REM  THE  5TH  SUBROUTINE  CALC 

ULATES  THE  ANGLE  THETA  WHICH 
IS  THE  ANGLE 

900  REM  FORMED  BETWEEN  THE  INTE 
RSECTING  PLANES  OF  THE  PLATE 
N’ S  TOP  AND  SLOPING 
910  REM  SIDES. 

920  GOSUB  6000 

930  REM  THE  6TH  SUBROUTINE  CALC 

ULATES  THE  THICKNESS  OF  THE 
[LATEN'S  BASE.  ITS 
940  REM  SQUARE  DIMENSIONS  WILL 
BE  IDENTICAL  TO  THE  DIAMETER 
OF  THE  SHAKER  HEAD. 

950  GOSUB  7000 

960  REM  THE  7TH  SUBROUTINE  CALC 
ULATES  THE  PLATEN  WEIGHT. 

970  GOSUB  8000 

980  REM  THE  8TH  SUBROUTINE  DISP 

LAYS  A  PLATEN  IN  HIGH  RESOLU 
TION  GRAPHICS 

990  REM  AND  LABELS  ALL  DIMENSIO 
NS  NECESSARY  FOR  FABRICATION 
EXCEPT  THE  TAPPED 
1000  REM  HOLE  LOCATIONS  FOR  TES 
T  SPECIMEN  MOUNTING. 

1010  GOSUB  9000 
1015  HOME 

1020  PRINT  "DO  YOU  WISH  TO  TRY  A 
NOTHER  DESIGN?" 

1030  PRINT  "IF  SO,  PRESS  KEY  <C> 

;  IF  NOT,  PRESS" 

1035  PRINT  "ANY  KEY.":  GET  C$ 

1040  IF  C$  *  "C"  OR  C$  *  "c"  THEN 

«5 
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1045  HOME 
1  047  VTAB  (15) 

1048  HTAB  (17) 

1049  PRINT  "GOODBY":  GOTO  2140 
2000  PRINT  "YOU  HAVE  ENTERED  A  R 

ESPONSE  THAT  THE" 

2010  PRINT  "COMPUTER  DOES  NOT  RE 
COGNIZE. 

2020  PRINT  "LETS  TRY  AGAIN." 

2030  PRINT 

2040  PRINT  "HAVE  YOU  COMPLETED  T 
HESE  PREL I M I  NARY" 

2050  PRINT  "DESIGN  DETAILS?  <Y  0 
R  N>" :  GET  AS 

2060  IF  A$  =  "Y"  OR  AS  =  "y"  THEN 
540 

2070  IF  AS  -  "N"  OR  AS  *  "n"  THEN 
21  00 

2080  IF  AS  <  >  "Y"  AND  AS  <  > 

"y"  AND  AS  <  >  "N"  AND  ASS  < 

>  "n"  THEN  21 00 

2100  PRINT  "PLEASE  PERFORM  THE  P 
RELIMINARY  LAYOUT" 

2110  PRINT  "STUOIES  REQUIRED  TO 
USE  THIS  PROGRAM." 

21 1 5  FOR  X  -  1  TO  2500:  NEXT  X 
2121  HOME  :  VTAB  (15):  HTAB  (17) 

:  PRINT  "GOODBY" 

2140  END 

2500  REM  THIS  SUB-ROUTINE  DETER 
MINES  WHICH  CORNER  HAS  THE 
GREATEST  SPECIMEN 
2510  REM  LOAD  BY  DIVIDING  THE  W 
EIGHT  IN  EACH  CORNER  BY  THE 
AREA  IT  DISPLACES. 

2520  HOME  :  PRINT 
2530  PRINT  "ASSIGN  EACH  CORNER  0 
F  THE  PLATEN  A" 

2540  PRINT  "NUMBER  SO  THE  WEIGHT 
S  OF  THE  AREAS  DIS-" 

2550  PRINT  "PLACED  BY  THE  TEST  S 
PECIMENS  AND  THEIR" 

2560  PRINT  "FIXTURING  CAN  BE  IDE 
NT  I F I  ED  BY  THEIR" 

2570  PRINT  "LOCATION  ASSIGNMENT. 

n 

2580  PRINT 

2590  FOR  X  *  1  TO  4 

2595  PRINT 

2600  PRINT  "WHAT  IS  THE  WEIGHT  0 
F  THE  TEST" 


Appendix  8 


41 


2610  PRINT  "SPECIMEN,  FIXTURING, 
BOLTS,  WASHERS, 

2620  PRINT  "ETC.  IN  CORNER  ";X;" 

?":  INPUT  W(X) 

2630  PRINT 

2640  PRINT  "WHAT  IS  THE  AREA  DIS 
PLACED  ON  THE" 

2650  PRINT  "PLATEN  BY  THE  TEST  S 
PECIMEN,  FIXTUR-" 

2660  PRINT  "ING,  ETC.  IN  CORNER 
" ; X ; " ? " :  INPUT  A(X) 

2670  PCX)  *  WCX)  /  A ( X ) 

2680  PRINT 

2685  PRINT  "THE  CORNER'S  LOAD  IS 
" ; P ( X ) ; "  PSI  " 

2690  NEXT  X 

2700  S  ■  PI :Y  -  1 

2705  IF  PC  1)  <  PC2)  THEN  PCI)  « 

PC 2) sY  »  2 

2710  IF  PCI )  <  PC3)  THEN  PCI  )  « 

P  (  3 )  :  Y  *  3 

27  20  IF  PCI)  <  PC  4)  THEN  PCI)  » 

PC  4 ) s  Y  *  4 

2730  PRINT  "CORNER  NUMBER  ";Yj" 
HAS  THE  GREATEST" 

2740  PRINT  "LOAD  OF  ";P(1);"  POU 
NDS  PER  SQUARE  INCH." 

2750  PRINT 

2770  PRINT  "PRESS  <ANY  KEY>  TO  C 
ONTINUE.":  GET  ZS 
2780  IF  ZS  -  >  CHRS  (0)  THEN  2 

790 

2790  HOME 
2800  RETURN 
3000  PRINT 

3010  REM  THIS  SUBROUTINE  IS  THE 
FIRST  ITERATION  TO  DETERM  I N 
E  WHAT  ANGLE  PHI 

3020  REM  SATIFIES  THE  EQUATION 
THAT  DETERMINES  THE  MINIMUM 
NATURAL  FREQUENCY. 

3030  C  -  CL  /  SQR  C2>)  -  (DIA  / 

2) 

3040  FOR  PHI  •  5  TO  45  STEP  .5 
3050  RAO  «  PHI  «  3.14159265  /  180 

3060  T  »  TAN  (RAD) 

3080  NU=  ( E  *  T  •  3 )  *  (  ( C  *  DEN 

*  T)  +  (12  *  PCI))) 

3090  OM  *  C  *  ((54  »  C  *  C  *  DEN  * 
DEN  *  T  »  T)  +  (8960  *  (PCI ) 

•  2))) 
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3100  IF  DM  <  *  0  THEN  3130 

3105  CF  *  48.47  *  ( ( NU  /  DM)  •  .5 
) 

3110  IF  CF  >  -  MF  THEN  3700 

3120  NEXT  PHI 
3125  SPEED-  1 50 

3130  PRINT  "THERE  IS  EITHER  AN  E 
RROR  IN  YOUR  INPUT" 

3140  PRINT  "DATA  OR  THE  MAXIMUM 
NATURAL  FREQUENCY" 

3150  PRINT  "OF  YOUR  PROGRAM  IS  T 
00  HIGH  FOR" 

3260  PRINT  "PRACTICALITY.  THE  HI 
GHEST  FREQUENCY" 

3270  PRINT  "POSSIBLE  FOR  THE  DES 
IGN  INFORMATION" 

3280  PRINT  "PROVIDED  IS  ";CF;"  C 
YCLES  PER  SECOND." 

3290  PRINT 

3300  PRINT  "IT  IS  SUGGESTED  THAT 
YOU  PERFORM  THESE" 

3310  PRINT  "CHECKS:" 

3320  PRINT 

3330  PRINT  TAB<  5)?"1.  RECHECK 
YOUR  OATA  INPUTS  TO" 

3340  PRINT  TAB<  8>;"SEE  IF  THEY 
ARE  CORRECT." 

3350  PRINT 

3360  PRINT  TAB (  5>;"2.  IF  THE  D 
ATA  IS  CORRECT  IN  #1" 

3370  PRINT  TAB  C  8) {"ABOVE,  THEN 
RECHECK  THE  LAYOUT" 

3380  PRINT  TAB (  8);"T0  DETERM  I N 
E  IF  THE  SPECIMENS" 

3390  PRINT  TAB (  8);"CAN  BE  SET 
CLOSER  TOGETHER.  THE" 

3400  PRINT  TAB (  8) {"NATURAL  FRE 
QUENCY  INCREASES  VERY" 

3410  PRINT  TAB (  8);"RAPIDLY  AS 
THE  LENGTH  OF  THE" 

3420  PRINT  TAB (  8>;"PLATEN*S  SI 
OES  DECREASES." 

3430  PRINT 

3440  PRINT  TAB (  5>;"3.  IF  NONE 
OF  THE  ABOVE  CAN  BE" 

3450  PRINT  TAB (  8 ); "CORRECTED, 
YOU  WILL  HAVE  TO" 

3460  PRINT  TAB<  8>;"ACCEPT  A  RE 
SONANT  CONDITION  IN" 

3465  PRINT  T AB (  8);"THE  PLATEN. 

•t 
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3468  SPEED*  255 
3470  PRINT 

3480  PRINT  "WHEN  YOU  PRESS  <AN Y 
KEY>" 

3485  PRINT  "THE  PROGRAM  WILL  TER 
MINATE." 

3490  GET  Z$ 

3600  IF  ZS  *  >  CHRS  (0)  THEN  1 

045 

3700  HOME 
3710  RETURN 

4000  REM  THIS  SUBROUTINE  ADDS  T 
HE  NECESSARY  ADDITIONAL  TH I C 
KNESS  FOR  TAPPING 
4010  REM  THREADS  INTO  THE  PLATE 
N  FOR  THE  BOLTING  ON  OF  TEST 
SPECIMENS.  IT  THEN 
4020  REM  CALCULATES  A  NEW  PLATE 
N  LENGTH  TO  FORM  WEDGE  EDGES 

4030  NT  *  C  *  T  +  TK 
4040  REM  THE  ABOVE  CALCULATION 
DETERMINES  THE  NEW  TOTAL  TH I 
CKNESS  OF  THE  PLATEN 
4050  REM  AT  DISTANCE  '  C*  ALONG 
THE  PLATENS  DIAGONAL. 

4060  REM  THE  NEXT  CALCULATION  D 
ETERMINES  A  NEW  "C"  BASED  UP 
ON  THIS  THICKNESS 
4070  REM  CALLED  "NC". 

4080  NC  *  NT  /  T 
4100  RETURN 

5000  REM  THIS  4TH  SUBROUTINE  RE 
CALCULATES  THE  ANGLE  PHI  BAS 
ED  UPON  "NC"  AND  "NT" 

5010  REM  EXCEPT  IT  IS  CALLED  PH 
12  TO  KEEP  IT  SEPARATE  FROM 
PHI.  THE  EQUATION  IS 
5020  REM  THEN  REARRANGED  AGAIN 
TO  UTILIZE  LINES  3050-3710  0 
F  THE  2ND  SUBROUTINE. 

5040  FOR  PH  1 2  *  PHI  TO  55  STEP  . 
1 

5050  RAD2  *  PH  1 2  *  3.1  41  59265  /  1 
80 

5060  T2  *  TAN  ( RAD2 ) 

5080  NU2  *  (E  *  T2  *  3)  *  C  C  NC  * 
DEN  *  T2 )  +  (12  »  P(1)>) 

5090  DM2  =  NC  *  ( ( 54  *  C  *  C  »  DE 
N  *  DEN  *  T2  *  T2 )  +  (8960  * 
(P( 1 )  •  2)  )  ) 
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5100  IF  DM2  <  *  0  THEN  3130 

5110  CF2  -  48.47  *  ( ( NU2  /  DM2)  • 
.5) 

5125  IF  CF2  >  ■  MF  THEN  5200 

5130  NEXT  PH  I  2 
5200  HOME 
5210  RETURN 

6000  REM  THIS  SUBROUTINE  CALCUl 
ATES  THE  ANGLE  ■THETA"  WHICH 
IS  THE  ANGLE 

6010  REM  FORMED  BETWEEN  THE  INT 
ERSECTING  PLANES  OF  THE  PLAT 
EN'S  TOP  AND  ITS 
6020  REM  SLOPING  SIDES. 

6030  BETA  =  45  *  3.1  41  5926  /  180 
6  03  5  NG  =  (NC  +  C  D I A  /  2))  *  SIN 
(BETA) 

6040  NT2  *  (NC  +  <  D I A  /  2))  *  T2 
6045  IF  NG  •  0  THEN  HOME  :  PRINT 
"A  DIVISION  BY  ZERO  ERROR,  R 
ECHECK  YOUR  WORK." 

6050  TETA  «  ATN  ( NT2  /  NG) 

6060  HOME 

6  07  0  RETURN 

7000  REM  THS  SUBROUTINE  CALCULA 
TES  THE  THICKNESS  OF  THE  PLA 
TENf S  BASE.  IT  OOES 
7010  REM  SO  BY  MULTIPLYING  THE 
ANGLE  THETA  BY  HALF  THE  D I  AM 
ETER  OF  THE  SHAKER 
7020  REM  HEAD.  IT  ALSO  CALCULAT 
ES  THE  TOTAL  THICKNESS  OF  TH 
E  PYRAMID  FROM  ITS 
7030  REM  BASE  (ITS  TOP  SINCE  IT 
IS  INVERTEO)  OOWN  TO  WHERE  I 
TIS  TRUNCATED. 

7040  BASEHEIT  ■  TAN  (TETA)  *  (D I 
A)  /  2 

7050  PYRAHEIT  -  TAN  (TETA)  *  L  / 

2  -  BASEHEIT 

7  06  0  HOME 
7070  RETURN 

8000  REM  THIS  SUBROUTINE  CALCUL 
ATES  THE  PLATEN’S  TOTAL  WEIG 
HT. 

8010  REM  IT  STARTS  WITH  THE  TOT 
AL  VOLUME  OF  THE  TOP  SQUARE- 
RETANGULAR  PARALLEL- 
8020  REM  PIPED. 

8030  TPVOL  ■  TK  «  L  •  2 
8040  REM  NEXT  IS  THE  CALCULATIO 
N  OF  THE  TOTAL  VOLUME  OF  THE 
PYRAMID  PORTION  OF 


8050  REM  THE  PLATEN 
8060  PI  -  ((  TAN  (TETA)  •  L  /  2)  « 
L  •  2)  /  3 

8063  P2  -  (BASEHEIT  *  D  •  2)  /  3 
8065  PVOL  -  PI  -  P2 
8070  REM  NEXT  IS  THE  CALCULATIO 
N  OF  THE  BASE  VOLUME 
8080  BVOL  -  BASEHEIT  •  0  •  2 
8090  REM  NEXT  IS  THE  TOTAL  VOLU 
ME 

8100  TTVOL  -  TPVOL  +  PVOL  +  BVOL 
8110  REM  FINALLY,  THE  TOTAL  PLA 
TEN  WEIGHT. 

8120  PLTWT  «  TTVOL  *  DEN 
8130  PYRAHEIT  *  (  INT  (PYRAHEIT  » 
100))  /  100 

8140  BASEHEIT  «  (  INT  (BASEHEIT  * 
100))  /  100 

8150  TETA  «  <  INT  ((TETA)  »  (180  / 
3.1  41  59))  *  100)  /  100 
8160  PLTWT  -  INT  (PLTWT) 

8200  RETURN 

9000  REM  THIS  SUBROUTINE  PRINTS 
A  DISPLAY  OF  THE  SIDE  AND  B 
OTTOM  V  IEWS  OF'  THE 
9010  REM  PLATEN  AND  PRINTS  OUT 
THE  FINAL  INFORMATION  REQUJR 
ED  TO  DESIGN  THE 
9020  REM  PLATEN. 

9040  HGR 
9045  HCOLOR*  3 

9050  PRINT  CHRS  <4)"BL0AD  GRAPH 

n 


10000  VTAB  (22) 

10010  PRINT  "LENGTHS:";  TAB (  11) 
;"L»";L;"  IN.";  TAB (  22);"D* 

" • D I  A • "  IN." 

10015  'PRINT  "TO  CONTINUE  PRINTIN 
G  DESIGN  INFORMATION" 

10016  PRINT  "PRESS  <ANY  KEY > . " :  GET 
CS 
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10017  IF  CJ  ■  >  CHRS  (0)  THEN 

1  0020 

10018  GOTO  10010 
10020  VTAB  (22) 

10025  PRINT  "TO  CONTINUE  PRINTIN 
G  DESIGN  INFORMATION" 

PRINT  "PRESS  <ANY  KEY  >  .  " 

IF  CS  *  >  CHRS  (0)  THEN 


1  0026 
1  0027 

1  0030 

10028  GOTO 


1  001  0 
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.  "v  *  -vi' 


■HUUUA 


10030  PRINT  "T*";TK;"  IN.";  T AB ( 

1  2) ;"P*";PYRAHEIT;"  IN.";  T  AB ( 
28) ;"B*";BASEHEIT;"  IN." 

10035  PRINT  "TO  CONTINUE  PRINTIN 
6  OESIGN  INFORMATION" 

10036  PRINT  "PRESS  <ANY  KEY>  . " 

10037  GET  CS 

10038  IF  CS  •  >  CHRS  (0)  THEN 

10040 

1  003  9  GOTO  10010 

10040  VTAB  (22) 

10050  PRINT  "ANGLE  THETA*" ;TETA; 

"  0 EG.";  TAB (  22); "PLATEN  NT 
.  ;  PLTWT;  "  LBS." 

10060  VTAB  (24) 

10260  PRINT  "WHEN  F I N I  SHED , PR E S S 
<ANY  KEY  >  .  " 

10270  GET  CS 

10280  IF  CS  -  >  CHRS  (0)  THEN 

1  0300 

10290  GOTO  10000 

1  0300  TEXT 

10310  RETURN 


